RSE 4P2/ERF 基因 家 族 鉴定 及 盐 Wy 迫 下 表 送 分 析 


EKA, mn) TWW JI? Hz 
(吉林 农业 大 学 林学 与 草 学 学 院 ， 长 春 130118) 
摘 Æ: AP2/ERF 转录 因子 在 植物 生长 发 育 和 响应 非 生物 胁迫 中 发 挥 着 重要 的 作用 。 为 探 
AL REESE (Aquilegia vulgaris) 中 AP2/ERF 对 盐 胁 迫 的 响应 ， 该 研究 基于 前 期 试验 获得 的 
盐 胁迫 下 转录 组 数据 ， 通 过 生物 信息 学 方法 筛选 欢 帮 斗 菜 4P2/ERF 家族 基 因 , 分 析 其 生化 
特征 、 保 守 基 序 、 系 统 进化 等 ， 并 对 其 在 盐 胁 人 迫 处 理 不 同时 间 的 根 与 叶 中 的 表达 量变 化 进行 
分 析 , 利用 qRT-PCR 技术 对 候选 基因 表达 量 进行 验证 ,结果 表明 :(1) 筛 选 出 86 个 4v4P2/ERF 
基因 ， 其 编码 的 蛋白 质 为 154~722 aa， 相 对 分 子 量 为 14 763.3~79 069.47 Da， 等 电 点 介 于 
4.49~9.68 之 间 ， 偏 酸性 蛋白 较 多 ， 均 为 亲 水 性 蛋白， 对 AVAP2/ERF EE EE 胞 定位 预 
测 ， 大 多 数 定位 于 细胞 核 。 O) 二 级 结构 以 无 规则 卷曲 和 o- 螺 旋 为 主 ， 均 具有 AP2 保守 结 
构 域 ， 有 两 个 高 度 保守 的 基 序 Motif 1 和 Motif 2。 (3) 在 盐 胁迫 下 ，71 个 AvAP2/ERF 基因 
表达 量 发 生变 化 ,其 中 叶片 中 差异 表达 基因 18 个 , 根 中 19 个 ; EOBEESE EU 7F AP2/ERF 
基因 聚 类 为 5 个 亚 家 族 ，15 个 亚 组 ， 通 过 表达 分 析 及 同 源 关 系 ， 确 定 3 个 响应 盐 胁 迫 的 基 
因 AvAP2/ERF-56, AvAP2/ERF-61 与 AvAP2/ERF-80, & qRT-PCR 结果 与 转录 组 数据 一 致 。 
ZH AKRA RRK K AP2/ERF 基因 的 功能 及 道 境 响应 机 制 莫 定 了 基础 。 
关键 词 : 欧 态 斗 菜 ，AP2/ERF， 生 物 信息 学 ， 盐 胁迫 转录 组 ， 表 达 分 析 
中 图 分 类 号 : Q943 文献 标识 码 : A 
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Identification of AP2/ERF gene family in Aguilegia vulgaris 


and expression patterns analysis under salt stress 


WANG Yufeng, MENG Yuan, YU Haihang, CUI Dingyuan, BAI Yun* 

(College of Forestry and Grassland Science, Jilin Agricultural University, Changchun 130118, China) 
Abstract: AP2/ERF transcription factors play important roles in plant growth, development and 
stress response. In order to explore the roles of Av4P2/ERF genes in salt tolerance, we identified 
AP2/ERF gene family in Aguilegia vulgaris transcriptome via bioinformatic methods, and then 
investigated their physical and chemical properties, conserved motifs, phylogenetic relations, and 
expression patterns of these genes in leaves and roots under salt stress. The expression of 
candidate genes was verified by qRT-PCR. The results were as follows: (1) 86 Av4P2/ERF genes 
were identified which encoded 154-722 amino acids, the molecular weight was 14 763.3-79 
069.47 Da and isoelectric point ranged from 4.49 to 9.68. Most of them were slightly acidic 
proteins and all of them were hydrophilic. Most of AvAP2/ERF were localized in nucleus. (2)The 
similarity of secondary structure was high, which was proportionally composed of random coil 
and a-helix. The members all contained AP2 domains, and two conserved motifs were predicted. 
(3) In different stages of salt treatment, 71 AvAP2/ERF genes responded to salt stress. There 
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werel8 and 19 differentially expressed genes in leaves and roots, respectively. 86 AP2/ERF genes 
of Aquilegia vulgaris were divided into 5 subfamilies clustering with Arabidopsis thaliana; 
AvAP2/ERF-56, AvAP2/ERF-61 and AvAP2/ERF-80 of them might be involved in salt resistance, 
and the qRT-PCR results were consistent with sequencing expression trends. This study provides a 
reliable reference for further research on the function and stress response mechanism of AP2/ERF 
gene in Aquilegia vulgaris. 

Keywords: Aguilegia vulgaris, AP2/ERF, bioinformatics, salt stress transcriptome, expression 


analysis 


AP2/ERF 是 广泛 存在 于 植物 中 的 一 类 转录 因子 超 家 族 ， 参 与 植物 的 生长 发 育 和 响应 外 
界 胁迫 。 该 家 族 成 员 均 含有 一 段 或 两 段 由 60—70 个 氨基 酸 组 成 的 AP2 保守 结构 域 ， 根 据 含 
有 AP2 结构 域 的 结构 或 数目 不 同 ，AP2/ERF 超 家 族 可 进一步 分 为 AP2、 ERF, DREB、 RAV 
及 Soloist 5 个 亚 家 族 〈 荀 艳丽 等 ，2020) 。 近 年 来 ， 许 多 植物 的 AP2/ERF 转录 因子 成 员 及 
功能 已 经 得 到 了 鉴定 与 验证 ， 拟 南 芥 (Arabidopsis thaliana) 中 鉴定 出 147 个 (Nakano et al., 
2006) , 水稲 (Oryza sativa) 中 163 个 (Akhter et al., 2011) , £RXZ (Vigna radiata) 中 186 
^^ (Chen et al., 2022) , 便 粒 小麦 (Triticum durum) 中 271 个 (Faraji et al., 2020) , A4% 
(Vigna radiata) 中 116 个 (Ran et al., 2022) 。 

AP2/ERF 转录 因子 在 植物 的 生长 发 育 、 应 对 外 界 各 种 生物 胁迫 、 非 生物 胁迫 中 都 发 挥 
着 重要 的 作用 (Feng et al., 2020) 。AP2 亚 家 族 主 要 参与 植物 的 生长 发 育 ，Jofuku 等 (1994) 
首次 从 拟 南 芥 中 分 离 出 AP2 转录 因子 ， 即 与 开花 生理 相关 。DREB 亚 家 族 和 ERF 亚 家 族 主 
要 在 植物 的 非 生物 逆境 胁迫 中 起 作用 《洪林 等 ，2020) 。RAV 亚 家 族 参 与 植物 响应 各 种 生 
物 和 非 生物 胁迫 过 程 〈Liu etal., 2021) 。 目 前 关于 Soloist 亚 家 族 的 研究 相对 较 少 ， 已 知 其 
参与 水 杨 酸 的 生物 积累 和 基础 防御 的 正 向 调节 王海波 等 ，2018) 。 

盐 胁迫 是 影响 植物 生长 发 育 的 最 重要 的 影响 因子 之 一 ，AP2/ERF 转录 因子 在 植物 应 对 
盐 胁 迫 过 程 中 ,发 挥 着 举足轻重 的 作用 。 拟 南 芥 41ERF1 的 表达 受 茉莉 酸 (JA) 、 乙 (ETO 
和 脱落 酸 CABA) 信和 号 相互 作用 的 控制 ， 作 为 信号 传递 的 中 枢 ， 在 拟 南 芥 的 耐 盐 调节 和 干 
旱 调 节 中 都 发 挥 着 重要 作用 (Cheng et al., 2013) 。 小 麦 TaERF-6-3A ìi SAW JW AWR Er R, 
抑制 抗 氧化 相关 基因 RD294 和 P5CS1 的 表达 来 降低 植物 的 耐 盐 性 ,起 到 了 负 调 控 的 作用 (Yu 
et al., 2022) 。 陆 地 棉 GhERF13.12 基因 受 ABA 信和 号 影响 ， 参 与 且 氨 酸 的 生物 合成 ， 在 拟 南 
芥 中 过 表达 , 能 增强 活性 氧 ROS) 清除 基因 的 表达 ,提高 拟 南 芥 对 盐 胁迫 耐 受 性 〈Lu et al., 
2021) . 早 柳 SmAP2-17 可 以 结合 SOS3 和 ABIS 的 启动 子 激活 它们 的 表达 , 在 盐 胁 迫 的 调控 
中 起 关键 作用 (Chen et al., 2022) 。 

KLE CAguilegia) 花卉 是 重要 的 宿根 花卉 ， 其 花形 独特 ， 花 色 艳 丽 ， 深 受 人 们 喜爱 
(M.Carmen et al., 2010; Jesus et al., 2015) ， 但 近年 来 ， 随 着 土地 盐 碱 化 的 加 剧 ， 炎 斗 菜 属 
花 齐 的 应 用 受到 限制 ， 盐 胁迫 使 其 观赏 效果 大 大 降低 。 课 题 组 前 期 试验 发 现 ， 欧 磷 斗 羔 
(Aquilegia vulgaris) 具有 很 强 的 抗 寒 性 和 耐 盐 性 ， 在 东北 地 区 可 以 露地 越冬 ， 管 理 方便 ， 
本 研究 以 欧 磋 斗 菜 盐 胁迫 处 理 不 同时 间 的 转录 组 测序 数据 为 基础 , ôl 29. 4P2/ERF 家族 基 因 , 
利用 生物 信息 学 方法 对 其 进行 分 析 ， 并 结合 试验 验证 ， 拟 探讨 以 下 问题 : (1) RE 
AP2/ERF 家 族 基因 编码 的 蛋白 质 的 理化 性 质 ; (2) 二 级 结构 及 保守 结构 域 特征 ; (3) 对 
盐 胁迫 不 同时 间 根 与 叶 中 4P2/ERF 基因 表达 模式 进行 分 析 、 系 统 进化 分 析 并 进行 分 类 及 功 
能 预测 ， 为 下 一 步 研究 欧 粮 斗 菜 4P2/ERF 基因 在 抗 盐 过 程 中 的 生物 学 功能 提供 参 


1 材料 与 方法 
1.1 转录 组 测序 试验 材料 


F 


N 
dk 
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期 试验 获得 , 种 子 为 吉林 农业 大 学 观赏 植物 资 
队 保 存 ， 于 2020 年 6 月 播种 ， 待 幼苗 长 至 6 片 真 叶 时 ， 选 择 长 势 良好 、 株 高 冠 幅 相 近 


的 植株 ， 进 行 200 mmol:LINaCl 溶液 处 理 ， 在 0h (CCK) 、12h、24 h. 48 h 后 分 别 对 根 和 


叶片 进行 取样 ， 重 复 3 次 ， 
后 置 于 -80 人 冰箱 保存 备 
分 析 。 


~ 
mi 


盐 溶液 浓度 及 处 理 时 间 均 


由 预 试验 确定 。 样 品 在 液 氮 中 迅速 冷冻 


1.20 KAX 4P2/ERF 基因 的 鉴定 及 生物 信息 学 分 析 


在 Tair Chttps://www.arabidopsis.org/) 上 捜索 


]， 由 北京 诺 禾 致 源 科技 股份 有 限 公 司 进行 转录 组 测序 并 进行 数据 


下 载 拟 南 芥 4P2/ERF 家族 的 基 因 序列 , 


作为 参考 序列 。 将 转录 组 序列 与 拟 南 芥 序列 进行 本 地 BLAST 比 対 , E 值 设置 为 1e-5， 将 得 


到 的 核酸 序列 通过 TBtools ( Chen et al, 2020 ) 翻译 为 蛋 


白质 序列 ， 


提交 NCBI 


(https://www.ncbi.nlm.nih.gov/) ~ SMART (http://smart.embl-heidelberg.de/ ) 预测 结构 域 ， 


保留 含 


Chttps://itol.embl.de/) 进行 美化 。 


完整 AP2/ERF 结构 域 的 蛋白 序列 。 利 月 
NJ) Mj ££ WX REGE 3E 5 11TH Zr APZ/ERF 基因 的 系统 发 育 
通过 在 线 软件 p/Mw (http://Web.expasy.Org/protparam/ ) 


进化 树 ， 


分 析 蛋 白质 的 分 子 量 、 氮 基 酸 数目 、 


A E 
于 


包 点 、 平 均 杀 水 性 等 理化 性 质 ， 利 月 


(http://cello.life.nctu.edu.tw/) 347 E 4 


胞 定位 。 利 月 


预测 和 蛋白质 
过 TBtools 进行 可 视 化 。 


二 级 结构 ，MEME (http://meme-suite.org/tools/meme ) 分 析 基 因 


1.3 WEBESESE APZ/ERF 基因 表达 模式 分 析 
将 AP2/ERF 家 族 基因 在 盐 胁 迫 0h (CK) 、12h、 
24h, 48h 的 表达 量 通 过 FPKM (fragments per kilobase million) 标准 化 ; AI 


IET OBESSE SER. Porc, 


H MEGA 7.0 软件 , 以 邻接 法 Cneighbor-joining， 


并 使 用 iTOL 


在 线 软件 MBC 
H Prabi Chttps://prabi.ibcp.fr/htm/site/web ) 
保守 基 序 ， 


并 通 


] DEGseq 软件 


进行 处 理 组 和 对 照 组 的 比较 , 差异 表达 显著 的 标准 为 Ilogz(Fold Change)|7 1. P(padj) 0.05, 


视 为 差异 基因 ( 牛 苏 燕 等 , 2022) 。 利用 TBtools 绘制 热 图 谱 , Row Scale 标准 化 , Cluster Rows 
聚 类 ， 其 他 参数 默认 ， 并 进行 表达 模式 分 析 。 
1.4 实时 荧光 定量 PCR 

为 了 验证 转录 组 测序 的 结果 及 目标 基因 的 表达 模式 , 设计 特异 引物 ( 表 1) , IPP2 作为 


内 参 基 因 (Sharma & Kramer, 2013) ， 经 实时 荧光 定量 PCR 检测 欧 炎 斗 菜 部 分 AP2/ERF 基 


因 在 盐 胁 迫 下 不 同时 间 在 根 和 叶 中 的 表达 。 月 


H SYBR Green I 检测 特异 引物 的 PCR 产物 ， 


反应 体系 如 下 : 2xSYBR Mix 主 混合 物 10 uL, EF 3514 1 uL, 模板 (cDNA) 2 uL, ddH2O 
6uL. PCR 反应 程序 : 预 变性 (95 'C, 1min) ， 放 大 定量 (95C, 20s; 60C, 20s; 72'C, 
30s, HE 40 次 ) ， 熔 解 曲线 (60 'C 295 CO 。 每 个 样品 3 个 独立 的 生物 学 重复 ，2 个 技 


术 重复 。 使 用 2G 法 定め 
或 下 调 ) 的 标准 化 比较 。 


基因 


的 相对 表达 量 ，log2(Fold Change) 进 行 基因 


表 1 欧 态 斗 菜 区 光 定 量 检测 引物 序列 


Table 1 
引 物 名 称 Primer name 
AvAP2/ERF-56F 
AvAP2/ERF-56R 
AvAP2/ERF-61F 
AvAP2/ERF-61R 
AvAP2/ERF-80F 
AvAP2/ERF-80R 
IPP2F 
IPP2R 


序列 ($-3 り Sequence (5'-3') 
GTATGGTGCCTCCCCTCGTT 
GCCCCTTGGTCTTGAACCTG 
GCGAAGTAGACGCAATGGACC 
AGCTGGCACTTTACGACGCT 
ATACGAAAGGCGGCAAGTGA 
CAACCCTGGCATTCCAAACTC 
CAGGTGAAGACGGACTGAAGTTATC 
CCAAGACTGGAAAAAAGACCACAC 


相对 表达 量 


(上 调 


i 


Sequence of the primers for fluorescent quantitative detection of Aquilegia vulgaris 


202301.00004v1 


chinaXiv 


2 结果 与 分 析 


2.1 转录 组 测序 结果 简介 
此 次 共 获 得 转录 组 文库 12 个 , 过 滤 后 每 个 文库 至 少 获得 41 521 502 个 clean reads、6.23 


G clean bases。 所 有 12 个 文库 中 的 Q30 均 超 过 94.81%, GC 均 在 41.51-42.72 之 间 ， 共 组 装 


ChinaXiv 合 作 期 刊 


序列 28 088 个 。 各 个 样本 比 对 到 基因 组 的 百分比 均 高 于 70%， 表 明 转 录 组 数据 与 参考 物种 


接近 。3 次 生物 学 重复 R2 均 大 于 0.85， 组 内 3 次 重复 存在 较 好 的 相关 性 。 为 了 进一步 挖掘 


盐 胁 迫 下 各 个 通路 的 基因 变化 ,以 P(padj) «0.05 为 筛选 条 件 , 筛选 出 13 个 显著 富 集 的 KEGG 


通路 ，2 197 个 差异 表达 基因 。 
2.2 WRJ 4P2/ERF 基因 家 族 鉴定 与 系统 发 育 分 析 
本 研究 从 欧 炎 斗 菜 转录 组 中 共 鉴 定 出 86 个 4P2/ERF 基因 家 族 成 员 ， 分 别 命 名 为 


AvAP2/ERF-1-AvAP2/ERF-86 . WW: Ej TU Fl JF 4P2/ERF 基因 构建 系统 发 育 进 化 树 ， 结 


果 如 图 1 所 示 ，Av4P2/ERF 分 布 在 ERF、DREB、AP2、RAYV 和 Soloist 这 5 个 亚 家 族 中 ， 
步 分 为 B11、B2、B3、B4、B5 fü B6 6 个 亚 组 ，DREB 亚 家 族 分 为 
A1、A2、A3、A4、AS 和 A6 6 个 亚 组 。 在 鉴定 的 86 个 4v4P2/ERF 中 , 15 个 属于 AP2 XE 
家 族 ， 占 总 数 的 17.44%; 37 个 属于 ERF 亚 家 族 ， 占 43.02%; 29 个 属于 DREB 亚 家 族 ， 占 
33.78%; 4 个 属于 RAV 亚 家 族 ， 占 4.6590; 还 有 1 个 属于 Soloist 亚 家 族 ， 占 1.16%。 


HF, ERF 亚 家 族 i 


WEE 86 个 4P2/ERF 基因 所 编码 的 蛋白 质 的 氨基 酸 数目 在 154~722 之 间 , 分 子 量 在 


ml) ummummmum 


1 Ml MTH TF 4P2/ERF 基因 的 系统 进化 树 
Fig.1 Phylogenetic trec of AP2/ERF genes in Aguilegia vulgaris and Arabidopsis thaliana 
2.3 欧 凌 斗 菜 AP2/ERF 蛋白 的 理化 性 质 分 析 与 亚 细 胞 定位 


14 763.3~79 069.47 Da 之 间 ， 等 电 点 介 于 4.49-9.68 之 间 ， 差 异 较 大 ， 暗 示 了 其 基因 功能 的 


多 欄 性 。 酸 性 


蛋白 数目 


(54) TRH 


蛋白 数目 


(32) , 平均 奈 水 性 的 値 費 小 干 0， 乡 为 亲 


水 蛋白 。5 个 定位 于 叶绿体 ，1 个 定位 于 线粒体 ，1 个 定位 于 细胞 质 ， 其 余 79 个 均 定位 于 细 


胞 核 ( 表 2) 。 定位 在 核 外 的 基 因 
AX CRITE, 2022) , Soloist | 
表 2 KYE AP2/ERF 的 理化 性 质 分 析 与 亚 台 


在 遗传 进化 上 
严 家 族 的 蛋白 
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EX 


Ua 


“rJ 


定 的 独立 必 


E， 与 家 族 内 其 他 基因 差异 


AvAP2/ERF-17 被 定位 在 细胞 质 。 


胞 定位 


Table 2 Analysis of physical and chemical properties and subcellular localization of AP2/ERF in 


AvAP2/ERF-1 

AvAP2/ERF-2 

AvAP2/ERF-3 

AvAP2/ERF-4 

AvAP2/ERF-5 

AvAP2/ERF-6 

AvAP2/ERF-7 

AvAP2/ERF-8 

AvAP2/ERF-9 

AvAP2/ERF-10 
AvAP2/ERF-11 
AvAP2/ERF-12 
AvAP2/ERF-13 
AvAP2/ERF-14 
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Table 3 Analysis of secondary structure of AP2/ERF in Aquilegia vulgaris 
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AvAP2/ERF-27 24.40 9.60 66.00 AvAP2/ERF-70 19.39 24.24 56.36 
AvAP2/ERF-28 34.12 12.17 53.71 AvAP2/ERF-71 19.48 14.94 65.58 
AvAP2/ERF-29 29.06 14.66 5628 AvAP2/ERF-72 19.08 16.76 64.16 
AvAP2/ERF-30 30.70 7.89 61.40 AvAP2/ERF-T3 24.58 11.3 64.12 
AvAP2/ERF-31 26.15 14.08 59.77 AvAP2/ERF-74 26.2 15.51 58.29 
AvAP2/ERF-32 25.08 9.15 65.76 AvAP2/ERF-15 23.35 7.00 69.65 
AvAP2/ERF-33 23.66 18.30 58.04 AvAP2/ERF-76 23.26 12.21 64.53 
AvAP2/ERF-34 18.06 12.74 69.20 AvAP2/ERF-TT 23.26 12.21 64.53 
AvAP2/ERF-35 19.67 10.39 69.94 AvAP2/ERF-78 30.74 15.58 53.68 
AvAP2/ERF-36 21.69 12.05 66.27 AvAP2/ERF-79 21.98 20.69 57.33 
AvAP2/ERF-37 25.00 11.53 63.47 AvAP2/ERF-80 32.47 16.02 51.52 
AvAP2/ERF-38 17.31 12.84 69.85 AvAP2Z/ERF-81 31.49 14.89 53.62 
AvAP2/ERF-39 43.95 12.11 43.95 AvAP2/ERF-82 25.43 14.45 60.12 
AvAP2/ERF-40 37.23 5.32 57.45 AvAP2/ERF-83 31.18 10.22 58.60 
AvAP2/ERF-41 35.90 10.26 53.85 AvAP2/ERF-84 29.68 13.24 57.08 
AvAP2/ERF-42 35.91 13.13 50.97 AvAP2/ERF-85 33.33 8.21 58.45 
AvAP2/ERF-43 36.65 14.91 48.45 AvAP2/ERF-86 21.33 20.00 58.67 


2.5 K4 AP2/ERF 蛋白 的 保守 基 序 分 析 

通过 MEME 进行 保守 基 序 分 析 ， 结 果 显 示 ，Motif 1 和 Motif 2 的 保守 程度 在 欧 磋 斗 菜 
AP2/ERF 中 是 最 高 的 ， 除 AvAP2/ERF-13 和 AvAP2/ERF-23 外 , 其 余 84 个 均 包 含 该 段 基 序 。 
而 其 他 类 型 结构 域 的 保守 程度 会 因 所 属 亚 家 族 的 不 同 而 有 差异 〈 图 2) 。 在 相同 亚 家 族 中 ， 

除 AP2 结构 域外 ， 还 存在 着 一 个 或 多 个 相对 保守 的 其 他 类 型 结构 域 ， 如 AP2 亚 家 族 有 两 段 
保守 基 序 Motif 4 和 Motif 13, RAV 亚 家 族 有 两 段 保守 基 序 Motif 9 和 Motif 15, 而 这 在 其 他 
[家 族 中 则 不 具备 (图 3) 。 
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Fig.2 Conserved motif of AP2/ERF in Aquilegia vulgaris 
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3 WOEE3ESIE AP2/ERF 基 序 分 布 


Fig.3 Distribution of motif of AP2/ERF in Aquilegia vulgaris 
2.6 WCRE APZ/ERF 基因 表达 模式 分 析 

KREE APO/ERF 基因 在 盐 胁迫 下 不 同时 间 在 根 和 叶 中 的 表达 数据 绘制 热 图 谱 并 聚 
类 分 析 , 将 FPKM 值 大 于 0.05 视 为 有 效 表达 , 并 且 其 值 越 大 说 明 表 达 水 平 越 高 (Moran et al., 
2011) ， 在 根 中 ， 有 效 表 达 基 因 69 个 ， 在 叶 中 ， 有 效 表达 基因 50 个 ， 且 表达 量 均 随 盐 胁迫 
时 间 增 加 呈现 一 定 的 变化 趋势 《图 4) 。 
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L、R 分 別 代 表 叶 片 和 根 部 ,」 0、12、24、48 为 处 理 时 间 (h)。 下 同 。 


L and R represent leaf and root, respectively; 0, 12, 24, 48 are the processing time (h). The same below. 
4 Rice APD/ERF 基因 盐 胁迫 下 表达 热 图 谱 
Fig.4 Heatmap of AP2/ERF genes expression in Aquilegia vulgaris under salt stress 


与 对 照相 比 ， 表 达 差 异 显著 的 标准 为 |log2(Fold Change))>1. P(padj)-0.05, 根 中 有 19 
个 差异 基因 ， 叶 中 有 18 个 差异 基因 ， 分 布 于 不 同 的 盐 胁迫 时 间 中 ，8 个 4P2/ERF 在 根 和 叶 
中 丝 为 差异 基因 。log2(Fold Change)>1 为 上 调 基 因 ，log2(Fold Change)<-1 为 下 调 基因 。 在 根 
中 ， 上 调 基因 有 15 个 ， 其 中 4v4P2/ERF-2-3-22-40-53-64 在 盐 处 理 12 h 时 表达 量 显 著 上 调 ， 
AvAP2/ERF-29-50-56-57-61-80 在 24 h 时 显著 上 调 ， 而 AvAP2/ERF-46-54-56-81 在 48 h 时 显 
x E. 下 调 基因 为 4v4P2/ERF-35-51-58-82， 显 著 下 调 的 时 间 也 不 尽 相 同 。 在 叶 中 ， 上 调 
基因 有 13 个 ，Av4P2/ERF-21-42-44-46-47-72-86 的 表达 量 在 12h 显著 上 调 ，AvAP2/ERF-80 
在 24h 时 显著 上 调 ，Av4P2/ERF-24-56-61-64-71 的 表达 量 在 48 h 显著 上 调 ， 
AvAP2/ERF-23-30-37-58-73 为 下 调 基 因 (图 5) 。 

根据 拟 南 芥 中 已 知 414P2/ERF、AP2/ERF 抗 逆 基 因 ( 表 4) o, BERRAR 
系统 进化 和 关系， 确定 AvAPZ/ERF-56-61-80 为 抗 盐 候选 基因 ， 且 这 3 个 基因 在 盐 胁 迫 处 理 后 
的 根 和 叶 中 ， 均 表现 出 显著 上 调 ， 可 见 其 积极 响应 盐 胁 迫 。 进 一 步 通 过 qRT-PCR 验证 上 述 
3 个 基因 在 盐 胁 迫 下 的 表达 模式 ， 由 图 6 可 知 , Av4AP2/ERF-56 在 叶 和 根 中 ， 均 随 着 盐 胁 迫 
时 间 的 增加 ， 表 达 量 逐渐 上 升 ， 在 盐 胁 迫 48 h 与 对 照 组 存在 着 显著 差异 ， 与 转录 组 测序 结 
果 基 本 一 致 , 4v4P2/ERF-61 在 叶 中 ， 随 着 盐 胁 迫 时 间 增 加 ， 表 达 量 逐渐 上 升 ， 在 48 h EU 
显著 , 在 根 中 ， 盐 胁迫 12. 24h 上 调 显著 ，48h 时 下 降 ， 与 转录 组 数据 一 致 ; AvAP2/ERF-80 
基因 在 叶片 中 ， 随 着 盐 胁迫 时 间 的 增加 ， 表 达 量 上 调 ， 在 24、48 h 时 上 调 显著 ， 在 根 中 ， 
盐 胁迫 12、24、48 h 时 均 检 测 到 表达 量 显著 上 调 ， 与 转录 组 测序 结果 一 致 。 
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图 5 WOBESESE AP2/ERF 差异 基因 盐 胁 迫 下 表达 热 图 谱 
Fig.5 Heatmap of AP2/ERF differential genes expression in Aquilegia vulgaris under salt stress 


表 4 4P2/ERF 家 族 同 源 基 因 在 拟 南 芥 中 的 功能 研究 
Table 4 Studies on the function of AP2/ERF family homologous genes in Arabidopsis thaliana 


WERE 
拟 南 芥 基 因 同 源 基 因 
Arabidopsis Homologous 
分 组 thaliana genes in 
Group genes Aquilegia 
vulgaris 
AP2 ANT AvAP2/ERF-37 
(AT4G377502 


Soloist APD1 
(ATAG13040 2 


AvAP2/ERF-17 


ERF-B3 ERF1 AvAP2/ERF-40 
(473C232402  AvAP2/ERF-42 
AvAP2/ERF-44 


AvAP2/ERF-47 


功能 


Function 


控制 根 细胞 数量 、 大 小 ， 负 调控 拟 南 芥 耐 盐 性 


Control the number and size of root cells; 


Negatively regulates salt tolerance in Arabidopsis 
thaliana 

在 水 杨 酸 防御 信号 途径 中 的 PADA 基 因 下流 
发 挥 功能 ， 正 调控 水 杨 酸 的 生物 合成 从 而 提高 
拟 南 芥 抗 逆 性 


It functions downstream of PAD4 gene in 


salicylic acid defense signaling pathway and 
positively regulates salicylic acid biosynthesis to 
improve stress resistance in Arabidopsis thaliana 
耐 基 、 耐 早 , JA. ET 和 ABA 信号 传递 中 枢 
Salt tolerance, drought tolerance, JA, ET and 


ABA signal transmission center 


1.50 
1.00 
0.50 
0.00 
-0.50 
-1.00 
-1.50 
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盐 胁迫 下 


表达 量 分 析 


Fig.6 Relative expression of partial AP2/ERF genes expression in Aquilegia vulgaris under salt 


3 讨论 与 结论 


本 研究 从 课题 组 前 期 获得 的 欧 磋 斗 菜 盐 胁 迫 转录 组 
化 特征 、 系 统 进化 及 盐 
AvAP2/ERF 可 被 分 为 15 个 4P2 基因 、29 个 DREB 基 
因 。 在 当前 研究 的 大 多 植物 中 ,ERF 亚 家 族 具 
严 家 族 成 员 最 少 (Nakano et al., 2006) 
严 家 族 成 员 在 数量 上 也 具有 相似 的 规 得 


基因 ， 并 对 其 家 族 成 员 到 


1 ^N Soloist 基 
AP2, RAV, 
种 的 AP2/ERF 


Soloist 】 


因 各 


Uu 


zET 


2022) ， 这 说 明 植 物 AP2/ERF 基因 
RILES AP2/ERF 转录 因 
Ran et al., 2022) , WEH f KRŠ- AP2/ERF 转录 因 


这 与 前 人 的 看 
构 较 为 复杂 ， 有 上 暗示 
说 明 
能 分 工 上 有 所 不 同 ， 在 
在 AP2/ERF 转录 因 


完结 果 一 臻 


LEZ 


Tuy. 


其 功能 的 多 样 | 
胞 质 中 合成 后 进入 细胞 核发 挥 作 月 
音 号 传递 的 不 同 阶段 发 挥 作用 。 
是 其 重 白 序列 的 


stress 


测序 数据 中 鉴定 了 86 个 AvAP2/ERF 


o 


H 


在 进化 上 有 共同 的 起 源 。 
子 家 族 成 员 的 分 子 量 、 等 


E. REK AP2/ERF 的 


H; 也 有 在 叶绿体 中 下 


胁迫 下 表达 模式 进行 了 分 析 。86 个 
因 、37 个 ERF 基因 


、4 个 RAV 基因 


及 


最 多 的 成 员 , 然 后 依次 是 DREB、 


欧 炎 斗 菜 与 拟 南 芥 、 水 稻 等 物 
(Akhter et al., 2011; 陈 悦 等 ， 


-存在 较 大 差异 ， 
子 家 族 的 结 
胞 定位 大 多 数 在 细胞 核 ， 
0 线粒体 中 的 ， 这 说 明 其 功 


电 点 和 杀 水 和 


y 


V2 


AH 


AP2 结构 域 的 高 度 保守 


EE 要 结构 特征 , 通 


け 


该 结构 域 包 含 60~70 全 気 基 酸 残 基 , 按照 3 个 B- 折 县 和 1 个 o- 螺 旋 方 式 形成 典型 的 三 维 结构 


( 洪 林 等 , 2020) . TEFKEEGLSE AP2/ERF 转录 因子 家 族 


高 的 ， 是 构成 AP2 结构 域 的 是 
常 与 转录 因子 
一 亚 家 族 基 因 
家族 有 两 段 
5 家族 
5 家族 


Y p 


同 亚 家 族 特异 的 保守 基 序 在 转录 调控 中 同样 发 挥 了 重要 上 


的 功能 相关 (Sakuma et al., 2002) 。 欧 
具有 相同 或 相似 的 保守 基 序 , 可 能 具有 相似 的 生物 学 功能 和 调控 途径 。 如 AP2 
^T di H* Motif 4 和 Motif 13, 构成 了 另 一 个 长 度 为 72 aa 的 AP2 结构 域 , RAV 


!, 


E 要 部 分 , XX WI ADR 


完 基 本 


FEJ AP2/ERF 转录 因 


Motif 1 和 Motif 2 是 保守 性 最 
相符 。 保 守 的 结构 域 与 基 序 通 
子 家族 中 , 同 


也 有 两 段 保 守 基 序 Motif 9 和 Motif 15, 是 B3 结构 域 的 习 
不 具备 ， 这 符合 AP2 亚 家 族 和 了 RAV il 


家 族 的 基本 特 生 


要 组 成 部 分 ,而 这 在 其 他 


E( 苟 艳丽 等 ，2020)〉。 不 


HEH. 


根据 与 已 知 功 能 基因 的 同 源 性 , 结合 表达 趋势 推测 目的 基因 的 功能 是 目前 研究 中 常用 的 
手段 ( 号 宇 辰 等 , 2022 ) 。 本 研究 中 将 欧 炎 斗 菜 86 个 、 拟 南 芥 147 个 4P2/ERF 基因 分 为 5 
个 亚 家 族 ， 各 家 族 成 员 功 能 类 似 。 通 过 同 源 比 对 发 现 ，AP2 亚 家 族 Av4P2/ERF-37 可 能 通过 
控制 根 细胞 数量 、 大 小 ， 负 调控 拟 南 芥 耐 盐 性 (Meng et al., 2015) o Soloist 亚 家 族 的 
Av4P2/ERF-17， 推 测 其 在 水 杨 酸 防 御 信 号 途径 中 的 P4D4 基因 下 游 发 挥 功能 ， 正 调控 水 杨 
酸 的 生物 合成 从 而 提高 植物 抗 道 性 〈Mrunmay et al., 2014) 。 预 测 ERF 亚 家 族 的 


AvAP2/ERF-40-42-44-47 在 抗 非 生物 胁迫 激素 信号 转 导 中 起 着 中 枢 作 月 


H CCheng et al., 2013 ) ; 


AvAP2/ERF-50 和 AvAP2/ERF-51 促进 抗坏血酸 合成 ， 提 高 植物 耐 盐 性 (Zhang et al., 2012) ; 
AvAP2/ERF-56 和 AvAP2/ERF-57 通过 ABN 介 导 的 ABA 信号 通路 ， 延 组 水 淹 诱 导 的 早衰 ， 
通过 表达 量 上 调 提 高 植物 的 耐 盐 性 和 耐 旱 性 (Liu et al., 2012) 。 DREB 亚 家 族 的 
AvAP2/ERF-79-80-81 可 能 参与 ROS 清除 ， 从 而 提升 植物 耐寒 性 与 耐 盐 性 (Kidokoro et al., 
20210 , AvAP2Z/ERF-61 通过 表达 上 调 提 升 植物 抗 盐 性 (Sowmya et al, 2011) , 
AvAP2/ERF-64-66-68 可 能 作用 于 光 和 乙烯 信号 通路 的 交汇 点 或 下 游 (Lin et al., 2008) , 
AvAP2/ERF-73-82-86 则 可 能 通过 降低 表达 量 ， 提 高 植物 抗 盐 性 〈Karaba et al., 2007) 。 

同一 基因 在 不 同 组 织 中 的 表达 可 能 存在 较 大 差异 (Zhou & Rajesh, 2021; 王丽娟 等 
2022) 。 本 研究 分 析 欧 炎 斗 菜 盐 胁迫 转录 组 表达 数据 发 现 ， 根 中 上 调 基 因 15 个 ， 下 调 3 
4 个 ; 叶片 中 ， 上 调 基 因 13 个 ， 下 调 基 因 5 个 ;在 根 和 叶片 中 省 上 调 的 基因 仪 5 S, EP 
调 的 仅 1 个 。 根 和 叶片 中 ， 差 异 基 因 有 明显 区 别 ， 表 达 差 异 显 著 检 测 到 的 时 间 也 不 尽 相同 
说 明 不 同 组 织 响应 盐 胁 迫 的 时 间 与 分 子 机 制 不 相同 ， 不 同 基 因 在 抗 盐 过 程 中 功能 也 不 同 。 
AvAP2/ERF-56. AvAP2/ERF-61 与 AvAP2/ERF-80 在 欧 态 斗 菜根 和 叶 中 均 显 车 上 调 ， 说 明 其 
积极 参与 了 抗 盐 过 程 ， 其 表达 量 随 盐 胁迫 时 间 发 生变 化 的 模式 有 所 不 同 , 也 体现 出 3 个 基因 
在 抗 盐 过 程 的 不 同 功 能 。 通 过 与 已 知 拟 南 芥 抗 盐 4P2/ERF 基因 同 源 比 对 发 现 ， 
AvAP2/ERF-56,. AvAP2/ERF-61 与 AvAP2/ERF-80 3 个 基因 极 有 可 能 在 盐 胁 迫 诱 导 下 提高 表 
达 量 ， 进 而 响应 及 抵抗 盐 胁 迫 。qRT-PCR 验证 发 现 三 个 基因 在 盐 胁迫 下 表达 模式 与 转录 组 
数据 一 致 。 

综 上 所 述 ， 本 研究 筛选 鉴定 86 NKEA AP2/ERF 基因 ， 对 家 族 成 员 进 行 了 详细 的 特 
征 分 析 与 进化 分 类 ， 结 合 转录 组 表达 数据 ， 对 4P2/ERF 基因 在 盐 胁迫 下 表达 模式 进行 了 总 
结 , 筛选 出 与 盐 胁 迫 相 关 的 候选 基因 3 个 , AvAP2/ERF-56. AvAP2/ERF-61 与 AvAP2/ERF-80, 
可 能 在 抵御 盐 胁 迫 中 发 挥 着 重要 作用 ， 但 具体 基因 功能 仍 需 开 展 后 续 试 验 进行 验证 。 
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